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We theoretically study three-dimensional topological semimetals (TSMs) with nodal lines pro¬ 
tected by crystalline symmetries. Compared with TSMs with point nodes, e.g., Weyl semimetals 
and Dirac semimetals, where the conduction and the valence bands touch at discrete points, in 
these new TSMs the two bands cross at closed lines in the Brillouin zone. We propose two new 
classes of symmetry protected nodal lines in the absence and in the presence of spin-orbital coupling 
(SOC), respectively. In the former, we discuss nodal lines that are protected by the combination 
of inversion symmetry and time-reversal symmetry; yet unlike any previously studied nodal lines in 
the same symmetry class, each nodal line has a Z 2 monopole charge and can only be created (anni¬ 
hilated) in pairs. In the second class, with SOC, we show that a nonsymmorphic symmetry (screw 
axis) protects a four-band crossing nodal line in systems having both inversion and time-reversal 
symmetries. 


The study of topological semimetals has recently 
drawn much attention from both the theoretical and 
the experimental communities. Topological semimetals 
(TSMs) are systems where the conduction and the va¬ 
lence bands have robust crossing points in fc-space, com¬ 
pared with normal metals where the two bands have a 
direct gap at each k. Compared with normal metals, 
the Fermi surface (FS) of an ideal TSM has a reduced 
dimension: in two dimensions, a normal metal has one¬ 
dimensional (ID) FS while a TSM has zero-dimensional 
(OD) FS; and in three dimensions, a normal metal has 
two-dimensional (2D) FS while a TSM has ID or OD 
FS. More importantly, the states near the FS are charac¬ 
terized by a nontrivial topological number. These unique 
features of FS in TSMs give rise to exotic properties such 
as the existence of Fermi arcs on the surfac^ and the 
chiral anomaly in the bullP^. In three dimensions, Weyl 
semimetalJi^and Dirac se mim etalhave been intensively 

studied both theoreticalljl^MlSl and experimentalljEll^. 

In Weyl semimetals, the two bands cross at an even num¬ 
ber of discrete points in the BZ, around which the bands 
are non-degenerate and disperse linearly in all three di¬ 
rections. In Dirac semimetals, both the conduction and 
the valence bands are twofold degenerate and cross each 
other at an odd or even number of points. Both systems 
belong to the class of topological nodal-point semimetals 
(TPSM). 

In three dimensions, there is another class of TSMs 
where the conduction and the valence bands cross each 
oth er at c losed lines instead of discrete points, i.e., nodal 
lineJUElHll! These nodal-line semimetals (TLSM) are in 
the midway between TPSMs and normal metals: (i) at 
exact half-filling, the FS is OD, ID and 2D in TPSMs, 
TLSMs and normal metals and (ii) the density of state 
scales as po oc (E — Ef)'^, po oc \E — Ef \ and a con¬ 
stant in TPSMs, TLSMs and normal metals, a fact from 
which one expects distinct electron correlation effects in 
the three classes. For example, the screening effect in 
these new metallic states have been discussecP^. Very 


recently, many theoretical proposals of materials for re¬ 
alizing TLSM have emerged, including the realization 
in graphene netwo rkj^, in CaaPj^, in LalNp^, and in 
Cu3(Pd,Zn)]N|2^*25|_ Jn all these w orks, and also in earlier 
theoretical model studiei^HsHlII^ the nodal line has the 
following properties: (i) unlike a Weyl node, a single line 
node can shrink to a point and vanish by continuously 
tuning the HamiltoniarP^l and (ii) its stability requires 
the absence of spin-orbital coupling (SOC), and upon 
turning on strong SOC, each nodal line is either split or 
gapped du e to t he hybridization between opposite spin 
component d^° l ^^ l It is natural to ask if there is another 
class of nodal lines with nontrivial monopole charges, and 
if there are nodal lines that are robust even in the pres¬ 
ence of SOC, i.e., four-band crossing lines. These open 
questions motivate us to develop a more comprehensive 
theory on TLSM. 

Our results are presented in two parts on TLSMs with¬ 
out and with SOC, respectively. In the first part, we re¬ 
visit systems with P and T in the absence of SOC, and 
find a new class of nodal lines that can only be created 
and annihilated in pairs, characterized by a new Z 2 topo¬ 
logical invariant. For a closed surface around the nodal 
line, we define a new Z 2 invariant protected hy P * T, 
classifying all nodal lines into two classes, with and with¬ 
out a Z 2 -charge, respectively. Nodal lines with Z 2 charge 
can only be created and annihilated in pairs, as the total 
charge of the BZ must be zero. Finite perturbation can 
make a nodal line with Z 2 charge shrink to an accidental 
nodal point, but cannot gap it. In the second part, we 
discuss systems with P, T and strong SOC. We show that 
if there is an additional twofold screw axis, a four-band 
crossing line, or a double nodal line (crossing between 
two doubly degenerate bands), can be protected on the 
boundary of the BZ. This is the first analytic proof of the 
symmetry protection of a four-band crossing line. We ap¬ 
ply the resultant theory to e xplain the double nodal line 
found in earlier model studied^Zlsoi on SrIrOa. 

The symmetry we consider is the composition of P and 
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T, OT P * T, an anti-unitary symmetry that preserves the 
momentum of a single particle. In a system without SOC, 
we have = 1, = 1 and [P,T] = 0, which imply 

(P * T)^ = -1-1. The action of P * T on the atomic or¬ 
bitals can hence be represented by complex conjugation 
{K) up to a basis choice. Therefore, the single particle 
Hamiltonian, P(k), is a real matrix at every k in BZ. 
A real, gapped Hamiltonian has Z 2 topological classifi¬ 
cation in both one and two dimensions, indicated by the 
first and the second homotopy groups of the projector 
onto the occupied band^ 



0{M + N) _ 0{M + N) 

0{M) ® 0{Ny 0{M) ® 0{N) 


( 1 ) 


where M and N are the numbers of the unoccupied and 
the occupied bands {M,N > 1). Since H{k) is real, at 
each k there is a real representation for all eigenstates 
of P(k). Therefore for each projector onto the occupied 
bands, P(k), 1 — 2P(k) is an 0{M + N) matrix that is 
invariant under any rotation within the occupied (unoc¬ 
cupied) space, i.e., an element of the quotient group. 

The Z 2 classification of P(k) in ID directly leads to 
protected nodal points in 2D and nodal lines in 3D, which 
have been studied in Refs. |26l32l34| . To see this, consider 
a closed path, i.e., a loop, in the 3D BZ, along which the 
Hamiltonian is gapped. The Z 2 invariant for the loop is 
simply the Berry’s phase for all occupied bands, quan¬ 
tized to either 0 or tt, corresponding to the trivial and 
the nontrivial classes respectively. If a loop belongs to 
the nontrivial Z 2 class, it cannot shrink to a point and 
vanish without crossing a singularity. In 3D BZ, this im¬ 
plies a line of singularities threading through the loop 
[see Fig0 a)]. Given H(k) continuous, each singularity 
is where the gap closes between the conduction and the 
valence bands, i.e., a nodal point, and a line of singular¬ 
ities is hence a nodal line. Given a nodal line in 3D, any 
loop that interlocks with the nodal line has Berry’s phase 
of TT, while the other loops have zero Berry’s phase. Thus 
we conclude that P*T topologically protects a nodal line 
in 3D in the absence of SOG. Unlike the monopole charge 
of Weyl nodes, the ID Z 2 invariant does not forbid a sin¬ 
gle nodal line from being annihilated or created locally 
in fc-space in 3D. Gonsider the following two band model 
as example 


H{k) = (m - k^aj; + k^a^. (2) 


FIG. 1: (a) A nontrivial Z 2 invariant (Berry’s phase of tt) of 
any loop in 3D BZ implies a nodal line (solid line) passing 
through the loop, (b) The 2D Z 2 invariant for a nodal line in 
3D BZ dehned on a sphere enclosing the line, (c) The evolu¬ 
tion of a nodal line with zero monopole charge as parameter 
changes in the model of Eq. (d) The evolution of a nodal 
line with nonzero monopole charge as parameter changes in 
the model of Eq. if 


in 3D. It implies that with P *T, H(k) on the surface of 
a sphere can be topologically nontrivial. In 3D fc-space, 
if i?(k) on the surface of a sphere belongs to the nontriv¬ 
ial class, the sphere cannot shrink to a point and vanish 
without meeting a singularity. Since P * T cannot sta¬ 
bilize a point node, the singularity is a nodal line inside 
the sphere. 

The 2D Z 2 invariant for nodal lines in 3D fc-space, 
defined on a surface enclosing the line [see Figj^b)], 
is a new topological invariant. It indicates the ab¬ 
sence/presence of the obstruction to finding a smooth 
and real gauge for the periodic part of the Bloch wave- 
functions (wavefunctions for short hereafter) of the oc¬ 
cupied bands. This invariant can be constructed based 
on this observation. Gonsider a sphere divided by an 
equator into two halves, say the northern and the south¬ 
ern hemispheres. The wavefunctions on each hemisphere, 
\uy {0, (j))) and |u;^(d, (j))), can be real and smooth because 
each hemisphere is contractible. At each point on the 
equator, denoted by the azimuthal angle (j), where the 
two hemispheres meet, the projectors to the occupied 
space must be equal 

n^occ. nGocc. 


H (k) has a nodal line on the kz = 0 plane of radius y/m 
if m > 0. As we change m from positive to negative, the 
nodal line shrinks to a point at the origin and vanishes 
[see c)]. The reverse process creates a single nodal 
line from the origin. This is a key difference between this 
nodal line and a point node in Weyl semimetal or some 
Dirac semimetal. In the latter cases, each point node has 
a monopole charge, and therefore can only be created or 
annihilated in pairs. 

The nontrivial Z 2 classification for real, gapped Hamil¬ 
tonians in 2D indicates a new classification for nodal lines 


because the Hamiltonian is smooth on the whole sphere. 
Therefore the matrix 

M™„((/.) = (n^(|,</>)|nf(|,</>)) (3) 

is an Nocc-ky-Nocc orthogonal matrix. For the orthogonal 
group, there is 

ni[OiNocc)] = Z 2 . (4) 

When M{(j)) is Z 2 nontrivial, there is an obstruction to 
defining a smooth gauge on the whole sphere. In Ap- 
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pendixj^ we show how to obtain this invariant without 
using a smooth gauge on hemispheres. 

The new invariant classifies nodal lines in 3D sys¬ 
tems with P *T into two classes: with and without a Z 2 
charge. A nodal line with a Z 2 charge can be considered 
a Z 2 monopole, which can only be created or annihilated 
in pairs. This is easy to prove by contradiction: if a 
Z 2 monopole is created (annihilated) locally in fc-space, 
the wavefunctions sufficiently away from this point are 
changed by a small amount, so the Z 2 invariant on a sur¬ 
face far away from the point is unchanged, contradicting 
the assumption that a Z 2 monopole is created (annihi¬ 
lated) within. 

For concreteness, we construct a 3D Hamiltonian 
around such a Z 2 monopole 


composed with a fractional lattice translation; the com¬ 
mutation relation involving a nonsymmorphic symmetry 
is generally ^-dependent due to the translation partPSHSl. 
At r, a nonsymmorphic symmetry can be treated as its 
point group component as long as the commutation re¬ 
lations are concerned; at BZ boundary, the fractional 
translation makes the group structure of nonsymmorphic 
symmetries different from any point group, and leads to 
new types of band crossings and high degeneracie^^. 

In this paper we show that the presence of P, T and a 
twofold screw axis protect double nodal lines (four-band 
crossing lines) on the BZ boundary in a 3D system with 
SOC. Unlike P that only acts in the real space, a twofold 
screw axis along z (Sz) acts in both the real space {x, y, z) 
and the spin space (sx,Sy,Sz) simultaneously: 


H{k) = QxSx + QyTySy + QzSz + mTxSx, ( 5 ) 

where and Si are Pauli matrices acting on two isospin 
degrees of freedom and q = k — ko is the momentum 
relative to the origin of the k-p expansion. The spectrum 
is given by 

P(k) = + {^qlPql±mY. (6) 

The band crossing can be found by solving P(q) = 0, 
yielding = 0 and q^ + qy = |to|, i.e., a nodal line on 

the a:y-plane of radius \/\m\. As m changes from positive 
to negative, the radius decreases and shrinks to zero at 
TO = 0 but increases again when to becomes negative [see 
Figj^d)]. An explicit calculation of the Z 2 charge of this 
nodal line is given in Appendix [B| 

We emphasize that the SU(2) rotation plays an impor¬ 
tant role in protecting a line node. When SU(2) is bro¬ 
ken, or there is SOC, the composition of inversion and 
time-reversal ensures double degeneracy at each k in BZ. 
Any crossing point between two doublet bands is hence a 
four-band crossing; yet it is easy to see that a four-band 
crossing is not protected in any 3D system without addi¬ 
tional symmetries. Without SU(2), the symmetry P *T 
satisfies (P * P)^ = — 1, because = — 1 in a spin half 
system . A generic four-band model with P * T in the 
presence of SOC is 

P(k) = flSx + hTySy + hSz + fiTxSy -|- f^TzSy, (7) 

where each is a function of {kx^ ky, kz), P*T — K{iTy). 
A band crossing requires five equations to be satisfied, 
namely, /i = /2 = /a = A = /s = 0, which is impossible 
in a 3D BZ without fine tuning. It is natural to ask if 
additional symmetries can protect a four-band crossing 
nodal line in 3D systems with SOC, and further if yes, 
what are they. 

Recent work shows that in the presence of nonsym¬ 
morphic symmetries, four-band crossings can appear at 
high symmetry OD nodal points on the BZ boundarjBi]. 
A nonsymmorphic symmetry is a point group symmetry 


^ , X / k-a -kb c, 

Sz : ix,y,z) ^ (-X + 2^ 

where /i, A = 0,1 denoting the shift of the axis, 
(/ia/4, A6/4, z), from the inversion center and a, b and 
c are the lengths of three basis vectors. Combining the 
twofold axis and inversion generates the another symme¬ 
try: 

Rz'. (a:,j/,z) ^ (x-y,?/-y,-z-f c/2) (9) 

(■^a:: -^z) ^ ^y: ^z)- 

P is a mirror plane (if /i = A = 0) or a glide plane (if 
/r or A is nonzero) located at z = c/4. Consider the 
commutation relation between P * T and Rz- In real 
space we have 

{x,y,z,t) (x-y,?/-y,-z-kc/2,t) (10) 

p*'!^ (-a;-f y,-y-k y,z-c/2,-t), 

{x,y,z,t) P*y {-x,-y,-z,-t) 

„ / MO A5 X 

(- 2 ;-y,- 2 /-y, 2 ; + c/2,-t), 

and in spin space 



^y j ^z ) 


( 5 “ 

'^yi ^z) 

(11) 



p*t: 

(^a:) ^y ; 

-Sz), 




p^t: 

( ^X-i ~ 

'^yi "^z) 





i:^x-> Sy, 

— Sz), 


which we 

find 





Rz * (-P * 

T) = 


-\b,c){P 

* P) * Rz 

(12) 


= g-'^k^+iyk^+iXky y * T) * P 


There are two planes (mirror invariant planes) defined by 
kz = 0 and fcj, = tt in the BZ that are invariant under 
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Rz, on which the commutation relations given by Eq.( 12) 
differ by a minus sign. 

On each mirror invariant plane, the bands can be la¬ 
beled by their respective Rz eigenvalues. In a system 
with SOC, we have 


Rl ■ V: z) ^ {x- y - Xb, z), 

^z') ^ ■) ^ z') ■ 


(13) 


or 


p2 _ _ rp _ _ ifj,kx-\-iXky 

— -^—fia, — Xb,0 — ^ 


(14) 




k ■ e 


The minus sign is because Rz is equivalent to a 7r-rotation 
along z in the spin space so R^ includes a 27r-rotation, 
giving a —1 for a spin-1/2 system. Therefore, each 
band at = 0 and kz = tt either has Rz eigenvalue 
5 + = or g_ = —g+. In the presence 

of SOC and P * T, bands are doubly degenerate, and 
the degenerate bands are related to each other hy P *T. 
Suppose at {k^, ky, kg), where fcp = 0, tt, a Bloch function 
|'0(k)) is an eigenstate of Rz with eigenvalue g+, then we 
consider its degenerate partner P * T|'(/'(k)) under Rz 

Rz{P*T)\'tP{k)) = 

= e-^'‘°g-P*T\'ijj(k)). 

At fco = 0, the degenerate bands have opposite Rz eigen¬ 
values, and two sets of such doublet bands generally anti¬ 
cross: the bands with the same Rz eigenvalue hybridize 
and avoid crossing [see Figj^c)]. At fcp = tt, however, 
the degenerate bands have the same Rz eigenvalue. In 
this case, two doublet bands with opposite Rz eigenval¬ 
ues may cross each other along a nodal line, making a 
symmetry protected four-band crossing line, or a double 
nodal line [see Figj^b)]. In Appendix]^ we revisit the 
3D fc.p-model in Fq. 0 in the presence of Rz in addition 
to P *T and show the presence of a double nodal line; 
and in Sec. IV, we write down a formal invariant for the 
double nodal lines. 

We apply the theory for double nodal lines to the case 
of the iridate SrIrOa. An eight-band tight-binding model 
consistent with all symmetries in the spac e group has 
been shown to exhibit the double line^^ESESl around point 
U in BZ, but a general symmetry analysis for arbitrary 
number of bands is missing to pin down which symme¬ 
tries of the little group at U are protecting the double 
nodal line, whereas the other symmetries may be bro¬ 
ken without opening a gap or splitting the double line 
into two single nodal lines and/or point nodes. The little 
group at U is generated by P plus two screw axes: 

P:{x,y,z) {-x,-y,-z), (15) 

Sy:{x,y,z) {-X + a/2,y+ b/2,-z + c/2), 

Sz-.{x,y,z) {-x,-y,z + c/2). 

Following our theory, we see that P, T and Sy {Sz) can 
protect double nodal lines on the ky = t: (kz = tt) plane. 


(c) 

FIG. 2: (a) A band crossing line on some mirror invariant 
plane at ^2 = 0 or kz = n, with an arbitrary cut along which 
the band structure is plotted, (b) The band structure along 
the cut in (a) on kz = tt, with the corresponding eigenvalues 
of Rz shown, where the degenerate states of a doublet band 
have the same eigenvalue, (c) The same band structure on 
kz = 0, where the degenerate states of a doublet band have 
opposite Rz eigenvalues, which can anti-cross with another 
doublet band due to the spin mixing enabled by SOC. 


Therefore, the double nodal line predicted in Ref. m on 
ky = TT plane is protected by P, T and Sy. The other 
symmetries, including Sz, S^ = Sy * Sz, Mz = Sz * P 
and Gb = P * Sx, may all be broken without gapping 
or splitting the double line. These statements may be 
tested in future experiments, as the crystalline symme¬ 
tries can be modified by applying epitaxial strain on thin- 
film sample^^. 

We briefly comment on the possibility of surface states 
in the TLSMs proposed. The protection of the nodal lines 
in this work requires the presence of P, as opposed to the 
Weyl semimetals where no symmetry (other than transla¬ 
tion) is required. Since an open surface always breaks the 
inversion, no protected surface states, strictly speaking, 
are associated with these TLSMs; however, as shown in 
Ref. |30I3^ ■ when the conduction and the valence bands 
are nearly symmetric, there is a nearly flat surface band 
bounded by the projection of the nodal lin^^. 

We propose two new classes of three-dimensional topo¬ 
logical semimetals with nodal lines in systems with and 
without spin-orbital coupling. Without spin-orbital cou¬ 
pling, or equivalently, with SU(2) spin rotation symme¬ 
try, inversion and time-reversal can protect a nodal line 
that carries a Z 2 monopole charge, independent of the 
previously known tt Berry’s phase. Nodal lines with a 
nonzero monopole charge can only be created or annihi¬ 
lated in pairs, while the nodal lines in previous studies 
can be singly created and annihilated. In the presence of 
SOC, we prove that inversion plus time-reversal are in¬ 
sufficient to protect any band crossing in 3D, and an ad¬ 
ditional nonsymmorphic symmetry (twofold screw axis) 
can protect a double nodal line, where two sets of doubly 
degenerate bands cross each other. We apply the theory 
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to SrIrOa, and identify symmetries required to protect 
the four-band crossing found in earlier model studies, and 
also symmetries that can be broken without gapping or 
splitting the double nodal line. 

* fangc@mit.edu 
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The surface fiat band proposed in SrIrOa is protected by 
the mirror symmetry Mz and an approximate chiral sym¬ 
metry specihc to this systerrP^. The chiral symmetry makes 
the conduction and the valence bands symmetric in energy. 


Appendix A: The explicit Z 2 invariant for a gapped 
and real Hamiltonian in 2D without using a smooth 
gauge 

In this Appendix, we explicitly construct the Z 2 invari¬ 
ant for a smooth mapping from a two-sphere to the man¬ 
ifold of all occupied states of a gapped and real Hamilto¬ 
nian, 0{M + N)/0{M) © 0{N), without using smooth 
gauges on the northern and the southern hemispheres. 

In this construction, we only need fix the gauges on the 
north pole and the south pole, where the wavefunctions 
of the occupied states are \uNP,i) and usp,i) respectively, 
where i labels the occupied bands and both wavefunc¬ 
tions are real. For each azimuthal angle (f) G [0,27r), 
there is a half longitude connecting the two poles, L(</>). 
Consider the following path-ordered integral 

W((/>) = Pexp[ f dP{9,(j))deP{0,(l>)d9], (1) 

Jo 

or, in the discretized form 

IF((?i)= lim Y[ P{jTr/N,(j)), (2) 

j=0,...,N 

where P{9,(j)) is the projector onto the occupied states, 
namely. 


P{Q,4')= X! \'^ii.^:4>)){Ui[9,4>)\. (3) 

i -1 , . . . ,iVocC 

The physical meaning of W{(f>) is a parallel transport be¬ 
tween the occupied space at the north pole and the oc¬ 
cupied space at the south pole. Therefore, the following 
matrix must be unitary 


= {usp,i\W{(j))\uNP,j)- (4) 


orthogonal transforms. Due to 7ri[0(A)] = ^ 2 , 

can either be trivial or nontrivial, and this Z 2 invariant 

of is the Z 2 invariant of H{9,(j)). 

Given a matrix function in 0{Nocc), how do we 

know if it is trivial or nontrivial? Consider a unit vec¬ 
tor in the real vector space of Ngcc dimensions: ng = 
(0,0,...,!)^. Acting M{(j)) on this vector results in a 
path in 


n{(j)) = M{(j))no. (5) 


Define 


^ n(0) - (n(((i) ■ no)no 
“ |n(0) - (n(((i) • no)nor 
a{(j)) = cos“^[no • 


( 6 ) 


Here n'(^) is the unit vector that is perpendicular to ng 
and inside the plane spanned by ng and and a{(j)) 

is the angle between n{(j)) and ng. Define the rotation 

R{(l),t) = cos[Q;((/))t] + sin[Q:(0)t]ngn'^((5i) (7) 


where t G [0,1]. Under these definitions, the orthogonal 
matrix 


0{(j),t) = R{(f),t)M{(l)), (8) 

is a smooth interpolation between M{(j)) at t = 0 and 
Mq{4>) at t = 1 that satisfies 

Mg(^)ng=no. (9) 

From Eq.([^, we see that Mo^cj)) is an orthogonal matrix 
that is block diagonal: the entries in last row and in the 
last column are all zero except the unity on the diagonal. 
Therefore, the first Aocc —1 columns and rows of Mg(^) is 
an orthogonal matrix in Nqcc ~ 1 dimensional real vector 
space. This block is defined as the new M{(j)) and we 
repeat the process. On each step, we smoothly connect 
M{(p) to Mg((()) that leaves the last unit vector invariant, 
then reduce the dimension of the matrix by one. This 
can proceed as long as dim{M) > 2, because the acting 
M{(j)) on ng gives a closed path in since 

7ri(S'"^^) = 0, this path can always shrink to a point 
smoothly. 

When is reduced to an 0(2) matrix, we define 

the U(l) number 

/I((^) = Mn(0)+iMi2(</)). (10) 

Due to the periodicity of M{<p), we can define the winding 
number of /3((^) 

2TTnu] =-i J (H) 


Due to the reality of {unp^sp) and of the Hamiltonian, The Z 2 invariant is the parity of rivv, or, 
is also real and hence orthogonal. Mij{(j)) is hence 
a smooth mapping from a closed path to the space of ^ — mod( Tlyj^ 2). 


( 12 ) 
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Appendix B: The Z 2 invariant for Eq.(5) 


The Zi invariant for any surface enclosing the origin 
can also be calculated for the model in Eq.(5). We will 
set m = 0 to calculate the invariant, and since it is a 
quantized, a small nonzero m does not change the re¬ 
sult. Given a sphere around the origin, the smooth wave- 
functions of the occupied bands on the northern and the 
southern hemispheres are 


10 

l«2^) 

kf) 


Q 6 6 

(cos sin - , — cos 2 ’ 2 

, , 0 . .6 6 rp 

(sm - sm 0, 0, — sm - cos 0, cos-) , 

/ & & , ^ 

(sm— cos - COS0,0, cos - sm0) , 

(0, cos - sm 0, — sm -, cos - cos 0) . 


The connection matrix on the equator is found to be 


M(0) = — COS0 — i sin 0(72. 


( 2 ) 


This path corresponds to the nontrivial element of 
7 ri[ 0 ( 2 )] as the winding number is - 1 - 1 , hence the Z 2 in¬ 
variant is nontrivial. 


Appendix C: Double nodal line in a four-band k ■ p 
model with P, T and a twofold screw axis 


Eqs.Q guarantee that /i _3 vanish on = 0 plane. But 
it takes fine tuning to make 02 , 4,5 vanish, i.e., to obtain a 
band crossing. This is consistent with our previous result 
that when [Rz,P*T] — 0 , there is no symmetry protected 
band crossing. If {R,P * T} = 0, or at k = (0,0 ,tt), 
Rz is represented by Rz = isx up to a gauge and the 
constraints on /i=i, 2 , 3 , 4,5 by Rz are 

fi{qx,qy,qz) = fi{qx,qy,-qz), ( 2 ) 

f2,3A,b{qx,qy,qz) = -/2,3,4,5(fe,9y,-<7z)- 


Eqs.® guarantee that 02 , 3 , 4,5 vanish on = 0 plane. 
Therraore, in order to obtain a band crossing on this 
plane, we only need fi{qx,qy,0) = 0 , the solution space 
to which is generically a ID line, i.e., a double nodal line. 
Appendix D: Topological invariant for a double 
nodal line 


According to the text, the double line crossing is 
formed by two doublet bands that have opposite eigen¬ 
values of R = P * Sz (where Sz can be replaced by any 
twofold screw axis). For a given double nodal ring on the 
kz = TT-plane, we can hnd two points inside and outside 
the ring on the same plane, denoted by A and B. At A 
(B), each band is doubly degenerate, where the two de¬ 
generate states have the same eigenvalue of R, labeled 
by ^n(A) We can hence define the following 

invariant 


We revisit the 3D fc.p-model in Eq.(7), in the pres¬ 
ence of Rz in addition to P *T. When the origin of the 
k.p expansion is at k = (0,0,0) (k = ( 0 , 0 , 7 r), we have 
[Rz,P*T] = 0 {{Rz,P*T} = 0). On At k = (0,0,0), 
the Rz is presented by Rz = iSy up to a gauge. Rz puts 
constraints on the forms of 0 i=i, 2 , 3 , 4 : 

fiA^x,qy.qz) = -fi,3{qx^qy^-qz)i ( 1 ) 

/2,4,5(93:5 — /2,4,5(9a:5 9y5 92)- 


c= n 

2nGocc 


^2n(A) 

6n(S)' 


( 1 ) 


Since R only has eigenvalues of 5 + and g- = —g+, C = il 
is a Z 2 invariant. When 0=1, the double nodal line 
between A and B is trivial and can be gapped, and if 
0 = — 1 , the line is protected. 



